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Abstract 

A  simple  capacitor  is  perceived  as  being  capable  of  discharge  (or  recharge)  at  high  rates,  limited  only  by  a  small  equivalent  series 
resistance.  However,  in  the  case  of  electrochemical  capacitors,  based  on  high  specific  area  porous  electrode  materials,  power  limitations 
arise  due  to  the  complex-distribution  of  electrolyte  internal  resistance,  coupled  with  double-layer  or  pseudo-capacitative  elements.  The 
present  paper  quantitatively  examines  both  numerically  and  experimentally  the  effects  of  internal  electrolyte  resistance  as  evaluated  by  the 
current-response  functions  for  porous  capacitor  electrodes,  generated  in  cyclic  voltammetry  (CV)  experiments  and  for  voltage  versus  time 
responses  in  dc  charge/discharge  regimes.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Electrochemical  capacitors  are  based  on  utilisation  of  the 
distributed  interface  of  conducting  materials  having  large 
specific  areas  such  as  carbon  felts,  cloths,  powders  and 
aerogels  [1],  and  also  conducting  metal  oxides  (RuOa) 
[2],  nitrides  [3],  or  conducting  polymers,  including  chemi¬ 
cally  or  electrochemically  generated  RuOi  deposits  on 
carbon  substrates  [4,5]. 

At  carbon  surfaces,  the  capacitance  arises  mainly  from 
electrostatic  double-layer  charging  [6]  while,  for  example,  at 
R11O2,  it  arises  mainly  from  redox  pseudo-capacitance  (2,7] 
coupled  with  a  smaller  double-layer  capacitance. 

It  is  important  to  emphasise  that  electrochemical  capa¬ 
citors  cannot  functionally  substitute  for  batteries;  rather  their 
applications  and  performance  behaviour  are  complementary 
to  those  of  rechargeable  batteries.  Elsewhere  [8],  we  have 
delineated,  in  some  detail,  the  similarities  and  differences 
between  batteries  and  electrochemical  capacitors  for  elec¬ 
trical  energy  storage  and  delivery.  Primarily,  electrochemi¬ 
cal  capacitors,  now  often  referred  to  as  “supercapacitors” 
[7],  provide  high-power  performance  but  at  low  energy- 
densities  and  also  with  the  advantage  of  large  multiple 
cycle-life  of  ca.  100,000  cycles  to  deep  depths  of  discharge, 
DOD.  On  the  other  hand,  batteries  provide  high  specific 
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energy-densities  but  at  relatively  low  operational  power- 
densities  and  with  limited  recyclabilities,  especially  at 
appreciable  DODs.  This  complementarity  of  performance 
has  led  to  various  proposed  hybrid  applications  employing 
tandem  configurations  of  battery  and  capacitor  modules 
enabling  optimisation  of  both  power  and  energy-density 
performance  [9]. 

The  relationship  between  operational  power-density  (P) 
and  simultaneously  achievable  energy-density  (E)  is  com¬ 
monly  represented  by  so-called  Ragone  plots  [10,11].  E 
tends  to  decrease  as  P  drain  is  increased.  In  a  capacitor,  this 
effect  is  mainly  due  to  internal  ohmic  (IR)  drop  which 
increases  with  current  (/)  drain  into  a  load,  lowering  the 
voltage.  In  a  battery,  it  is  due  to  kinetic  polarisation  effects 
in  addition  to  IR  drop,  together  with  significant  effects 
arising  from  concentration  gradients  developed  on  charge 
or  discharge.  In  operation  of  a  capacitor  device,  for  funda¬ 
mental  reasons,  the  Ragone  relations  are  importantly  depen¬ 
dent  on  state-of-charge  [11]  (SOC)  since  the  operational 
voltage  ( V),  is  directly  proportional  to  charge  (Q),  V  =  QIC, 
so  that  both  P  and  E  are  determined  by  SOC.  In  battery 
performance,  related  effects  can  arise  but  are  secondary  and 
smaller. 

In  electrochemical  capacitors  employing  non-aqueous 
electrolytes,  the  internal  electrolyte  resistance  is  usually 
appreciably  higher  than  for  aqueous  solutions  but  this  is 
compensated  by  higher  operational  voltage,  which  enters 
“squared”  into  the  energy  density. 
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A  pure  capacitance,  such  as  that  of  the  double-layer  [6], 
can  be  discharged  at  high  rates,  limited  only  by  a  (usually 
small)  equivalent  series  resistance  (ESR).  However,  in  the 
case  of  electrochemical  capacitors  based  on  high-area, 
porous  materials,  an  unavoidable  distributed  resistance 
(R)  arises  (mainly  of  the  electrolyte  in  the  pores)  that  is 
coupled  in  a  complex  way  with  the  distributed  interfacial 
capacitance  elements  (C),  giving  rise  to  a  distribution  of 
“RC”  time-constants  and  a  resultant  power  spectrum.  This 
situation  was  treated  in  mathematical  detail  by  de  Levie 
[12,13]  in  a  series  of  papers.  The  important  general  con¬ 
clusions,  relevant  to  porous  capacitor  behaviour,  were:  (a) 
that  a  so-called  “penetration  effect”  (into  the  porous  matrix) 
applies  to  any  time-dependent  electrical  signal  —  current 
pulse,  potential  pulse  or  amplitude  voltage  (a.v.)  modulation 
—  so  the  charging  or  discharging  response  function  is  time- 
dependent;  and  (b)  that  the  overall  behaviour  under  a.v. 
modulation  is  a  transfer  function  characterised  by  a  —45° 
phase  angle  of  the  impedance  vector,  represented  in  the 
complex-plane. 

These  aspects,  treated  by  de  Levie  [12]  and  later  by  Keiser 
et  al.  [14]  for  pores  of  various  shapes,  provide  the  founda¬ 
tions  for  consideration  of  power  limitations  in  the  operation 
of  electrochemical  capacitor  devices  based  on  high-area, 
porous  electrode  materials. 

Because  of  the  relationship  between  the  response  cur¬ 
rents,  /,  to  repetitive  linear  voltage  modulation  (cyclic 
voltammetry  (CV))  at  sweep  rates  s  =  ±dV/dr,  giving 
directly  the  capacitance  C  =  I/s,  CV  has  become  a  major 
procedure  for  evaluating  performance  of  capacitor  devices. 
It  is  convenient  for  determining  cycle  life,  and  from  analysis 
of  the  shapes  of  the  response-current  voltammograms  as  a 
function  of  s,  it  provides  information  on  internal  resistance 
effects  and  consequent  dissipative  losses. 

The  present  paper  examines  consequences  of  such  effects 
in  quantitative  ways,  through  simulation  calculations,  espe¬ 
cially  for  cases  where  the  performance  characteristics  are 
evaluated  by  means  of  CV  and  de  charging  curves.  The 
results  of  simulation  calculations  for  various  conditions 
applied  to  voltammetry  and  dc  charging  behaviour  at  porous 
capacitor  electrodes  are  compared  with  experimental  results 
obtained  at  specially  fabricated,  non-aqueous  electrolyte 
capacitor  modules  containing  sub-optimal  electrolyte  con¬ 
centrations  in  order  to  enhance  internal  resistance  effects, 
thus,  enabling  them  to  be  more  quantitatively  characteriz- 
able. 


2.  Treatment  of  cyclic  voltammetry  (CV)  and  dc 
charging  at  porous  electrodes 

Equations  have  been  developed  for  potentiostatic  and 
galvanostatic  charging  transients  at  porous  and  tubular 
electrodes  of  finite  length  in  the  absence  of  appreciable 
Faradaic  currents,  for  systems  of  constant  electrolyte  resis¬ 
tance  [15,16].  The  case  of  linear  or  triangular  voltage  sweep 


(CV )  in  the  presence  of  significant,  but  constant,  IR  drop  for 
porous  electrodes  was  considered  in  detail  by  Austin  and 
Gagnon  [17], 

Effects  of  internal  electrolyte  resistance  and  its  changes 
under  certain  conditions  during  charging  and  discharging  of 
porous  capacitor  electrodes  are  examined  here  for  cases  of 
charge/discharge  under  CV  and  for  dc  constant  current 
cycling.  The  required  mathematical  equations  are  developed 
and  their  solutions  are  evaluated  with  numerical  simulations 
for  various  conditions. 

The  assumptions  made  in  the  calculations  include  the 
following: 

1 .  that  the  double-layer  capacitance,  C,  is  constant  over  the 
voltage  window  of  operation; 

2.  only  capacitative  processes  occur,  i.e.  there  is  no 
Faradaic  charging; 

3.  only  ohmic  losses  are  considered; 

4.  the  electrolyte  conductivity  is  either  (1)  constant,  or  (2) 
a  linear  function  of  the  local  potential  at  the  working 
electrode. 


3.  Governing  equations 

The  current  response,  I(Vm)  at  a  potential  lj„,  for  a  purely 
capacitative  process,  involving  a  constant  capacitance,  C,  is 
given  by: 

/(V„)  =  c(§)  (1, 

where  Vm  is  the  measured  potential  and  Vc  the  local  potential 
at  the  capacitative  electrode,  Ve  differs  from  Vm  by  the  “IR” 
loss  due  to  potential  dependent  passage  of  response  current 
when  the  system  includes  a  significant  series  resistance. 

3.1.  Constant  electrolyte  resistance 

The  local  electrode  potential,  during  the  anodic  direction 
of  polarisation  of  a  CV  experiment,  between  the  applied 
potential  limits  of  and  V2,  assuming  the  presence  of  a 
constant  electrolyte  series  resistance,  RQ,  is 

Ve  =  Vi+st-  IR0  (2a) 

and  during  the  cathodic  direction  of  the  CV,  between  the 
applied  potential  limits  of  V2  and  V\,  is 

Vc  =  V2  -  st  -  IR0  (2b) 

where  ,v  is  the  applied  potential  sweep  rate,  dV/d t.  Differ¬ 
entiating  Eqs.  (2a)  and  (2b)  with  respect  to  t,  and  substituting 
for  dV/dr  in  Eq.  (1),  results  in  the  following  first-order 
differential  equation  for  anodic  current  /a(Vm) 

4(Vm)  =  Cs-RqC—  for  0  <  t  <  (3a) 

at  s 
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and  for  cathodic  current  /c(Ve) 

/c(Vm)  =  -Cy-fl0C^ 

(V2-Vi)  2(V2-V1) 

for  — - -  <  t  <  — - -  C3b) 

s  s 

The  analytical  solutions  of  Eqs.  (3a)  and  (3b),  subject  to  the 
boundary  conditions 

,4.) 

and 


h  =  sC  — 


2  sC 


for  0  <  t  < 

and 

Ic  =  sC 


1+exp  ((Vi-V2)/sRoC) 
(V2  -  Vi) 


exp 


2  sC 


x  exp 


1  +  exp((y1  -  V2)/sR0C) 
(V2  —  V\)/s—  t 


RoC 


for  (V>-VQ<t<2<y>-V'). 


R0c) 


(5a) 


(5b) 


3.2.  Case  for  electrolyte  being  consumed  during  charging 
of  the  capacitor 

Consider,  next,  the  case  in  which,  during  the  course  of 
charging  of  the  capacitative  electrode,  the  electrolyte  is 
“consumed"  and,  as  an  approximation,  the  electrolyte 
resistance  then  becomes  a  linear  function  of  the  local 
electrode  potential,  thus 

R  =  R0  +  rVe  (6) 

In  real,  two-electrode,  high-area  capacitor  device  this  “con¬ 
sumption”  of  electrolyte  from  the  bulk  solution  is  a  measur¬ 
able  effect  and  arises  because  of  differential  cation/anion 
adsorption  between  the  two  electrodes  as  the  their  double¬ 
layers  become  charged  [18,19], 

As  a  result  of  Eq.  (6),  the  local  potential,  Ve,  during  the 
anodic  direction  of  polarisation  in  a  CV  experiment,  con¬ 
ducted  between  the  applied  potential  limits  of  Vm  =  14  and 
Vm  =  V2,  is 

yt  =  Vl+ st-  IR0  -  IrVe 

or  14  =  (Vj  +  st  -  IR0)(l  +  Ir)~l  (7a) 

while  during  the  cathodic  direction  of  a  CV  experiment, 
between  the  applied  potential  limits  of  V2  and  14,  Vc  is 

14  =  (214  -Vi  -st-  IRo)  ( 1  +Ir)~l  (7b) 


After  differentiating  Eqs.  (7a)  and  (7b)  with  respect  to  t  and 
substituting  into  Eq.  (1),  the  following  expressions  for 
current  as  a  function  of  time  result: 


/»  = 


C?(l  +  Ir)  —  C(Rq  +  rVi  +  rst)(dl/dt) 


for  0  <  t  < 


(1  +  Ir)2 
V2-Vi 


and 


(8a) 


— Cj(1  +  Ir)  —  C(Ro  +  2r!4  —  rV\  —  rst)(dl/dt ) 

(1  +Irf 


for 


Vi-Vi  <t<  2(14  -  Vi) 


(8b) 


s  s 

A  numerical  solution,  subject  to  the  boundary  conditions 
of  Eq.  (4a)  and  (4b),  may  be  determined  using  the  Euler- 
Cauchy  method  [20],  since  the  change  in  current  with  time 
may  be  expressed  as 


d 7a  _  Cslr  +  Cs  —  I  —  2 f~r  -  /  V2 
dt  CRq  +  CrV i  +  Crst 

Vo  -  14 

for  0  <  r  <  — - -  (9a) 

s 

and 


d  Ic  — Cslr  —  Cs  —  2 12r  —  /  V2 

dr  CRq  +  2CrV2  —  CrV  i  —  Crst 


v^-vi<t<2(y^-v}i 


(9b) 


with  an  initial  value  of  /a(r  =  0)  being  assumed.  The  numer¬ 
ical  solution  determined  using  Eqs.  (9a)  and  (9b),  for  given 
Rq  and  r  =  0  was  identical  to  the  analytical  solution  calcu¬ 
lated  using  Eqs.  (5a)  and  (5b). 


3.3.  Case  for  electrolyte  being  consumed  during  constant 
current  charging  of  the  capacitor 


The  electrical  response  during  charge/discharge  is  gov¬ 
erned  by  Eq.  (1),  but  becomes  simplified  in  the  case  of 
constant  current  to: 

AVe  =  ^At  (10) 

Given  that  the  measured  voltage  includes  the  IR  component, 
it  may  be  calculated  simply  as 

Em  =  14(1  +lr)  +  /R0  (11) 

i.e.  Vjn  is  determined  both  by  the  intrinsic  resistance,  R0,  of 
the  electrolyte  medium  and  its  change  with  polarisation 
potential  expressed  by  r  in  Eq.  (6). 


4.  Experimental 

4.1.  Comparisons  with  behaviour  of  real  capacitor  modules 

In  order  to  examine  the  validity  of  the  simulation  treat¬ 
ments  developed  in  Section  3,  comparisons  were  made 
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experimentally  with  the  charge/discharge  behaviour  of  two 
specially  fabricated,  non-aqueous  porous  carbon  electrodes, 
described  below.  The  experimental  test  procedures 
employed  CV  conducted  at  various  sweep-rates  and  dc 
charge/discharge  cycling  at  various  currents. 

The  temperature  for  all  experiments  was  25  ±  1°C.  Elec¬ 
trical  connections  to  the  modules  from  the  potentiostat  were 
made  in  a  two-electrode  configuration,  with  one  of  the  two 
capacitor  electrodes  serving  both  as  a  counter  and  reference 
(positive  terminal). 

4.2.  Capacitor  module  fabrication 


Type  14-50  (size  R6/AA)  hardware  was  used  to  construct 
the  sealed  single-cell  electrochemical  capacitor  modules 
used  in  these  tests.  The  porous  electrode  modules  employed 
1.0  and  0.40  M  tetraethylammonium  tetrafluoroborate 
(TEATFB)  in  dry  propylene  carbonate  as  electrolyte,  as 
described  elsewhere  [18], 

5.  Results  and  their  interpretations 

5.1.  Calculated  CV  responses 

5.1.1.  Case  of  constant  resistance  (Ro) 

In  Fig.  1,  current  responses  are  conveniently  plotted  as 
reduced  currents  (l(Vm)/sC)  as  a  function  of  measured 
potential  (Vc  +  IRq),  of  a  capacitative  electrode  exposed 
to  electrolytes  having  various  constant  electrolyte  resis¬ 
tances,  Rq,  between  0  and  50  Q,  (C=  1  F,  s  =  20  mV  s-1, 
Vi  =  0  V,  V2  =  1  V)  as  previously  discussed  in  some  detail 
[21],  For  sufficiently  small  resistance,  R0,  the  current 
response  is  essentially  that  of  a  pure  capacitor  (Fig.  1  for 
Rq  — >  1  Q).  Increasing  R0  increases  the  IR0  component  and, 
given  a  sufficiently  large  R0,  the  overall  behaviour  of  the 
resistor/capacitor  network  becomes  dominated  by  the  resis¬ 
tive  component.  Note  also,  that  the  form  of  the  CV  calcu¬ 
lated  for  the  constant  resistivity  simulation  is  symmetrical 
with  respect  to  an  in-plane  rotation  of  1 80°  which  renders  the 
IJsC -1  versus  V  response  superimposable  with  the  reverse 
scan.  This  indicates  that  the  beginning  of  discharge  behaves 
identically  to  the  beginning  of  charge  with  the  exception  of 
reversal  in  the  direction  of  current. 

Examination  of  the  reduced  current  as  a  function  of  the 
local  electrode  potential,  Ve,  (Fig.  lb)  illustrates,  given  a 
sufficiently  large  resistance  and  resulting  IR0  loss  compo¬ 
nent,  how  the  applied  and  actual  potentials  are  significantly 
different  so  that  the  actual  sweep  rate  at  the  working 
electrode  and  the  effective  range  of  applied  potential  is 
far  from  what  the  experimenter  intended.  As  well,  the  local 
sweep  rate,  dVe/dr,  for  large  R0,  is  not  constant,  and  the 
I/sC~l  versus  Vc  response  is  not  of  the  expected  rectangular 
form  for  a  purely  capacitative  device. 

Correspondingly,  similar  behaviour  results  with  increas¬ 
ing  capacitance,  C,  and/or  sweep  rate,  s,  while  maintaining  a 


Fig.  1.  Calculated  cyclic  voltammograms  expressed  as  reduced  response 
current  ( ItsC)  as  a  function  of  (a)  measured  potential  (Ve  +  IRo)  and  (b) 
local  working  electrode  potential  (Ve)  calculated  using  Eqs.  (5a)  and  (5b) 
for  C  =  1  F,  Vi  =  0  V,  V2  =  I  V,  s  =  20  mV  s1.  (if)  Indicates  direction  of 
increasing  R0  (1,  5,  10.  25  and  50  Q);  (  >  )  direction  of  voltammetric 
sweep. 


given  fixed  resistance  value.  In  the  former  case,  an  increase 
in  device  capacitance,  C,  results  in  an  increase  in  the  current 
/passed  (Eq.  (1)),  and  hence  an  increase  in  the  IRQ  losses.  In 
the  latter  case,  at  low  sweep  rate,  I  is  small,  so  that  IRq  losses 
are  negligible,  but  as  the  applied  sweep  rate  is  increased,  the 
current  increases,  so  that  IR0  losses  become  significant.  The 
dependence  of  “effective”  capacitance  on  sweep  rate  is 
of  particular  importance  as  it  provides  an  indication  of  the 
rate  performance  of  the  electrode,  i.e.  a  high  sweep  rate 
corresponds  to  a  high  rate  of  charge  or  discharge  of  the 
capacitor  electrode,  corresponding  to  greater  power  levels 
of  operation. 

Importantly,  these  results  show  that  even  though  the  actual 
capacitance  of  the  working  electrode  is  constant  and  inde¬ 
pendent  of  either  Ro  or  s,  the  “effective”  capacitance,  which 
is  determined  by  the  IR0  component,  is  not  constant.  Rather, 
the  effective  capacitance  decreases  with  increasing  R0  and/ 
or  s,  and  is  typically  much  less  than  C.  Thus,  less  and  less  of 
the  “theoretical”  charge  capacity  is  accessible  with  increas¬ 
ing  Rq  and/or  s,  leading  to  diminution  of  energy-density 
and  power-density  of  the  device  modelled  by  the  above 
equations. 
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5.1.2.  Case  where  R  a  linear  function  of  the  local  working 
electrode  potential 

We  now  examine,  in  some  detail,  the  results  of  simulation 
calculations  for  a  capacitative  electrode  subject  to  increasing 
electrolyte  resistivity  with  increasing  state-of-charge,  i.e. 
subject  to  a  resistance  component  dependent  on  Ve  as 
represented  by  Eq.  (6).  The  effects  of  the  variable  resistance 
on  the  CV  response  for  three  distinct  cases  are  examined:  (a) 
constant  s,  C  and  r  (the  voltage-dependent  component  of 
resistance),  and  various  values  of  the  so-called  “fixed” 
component  of  the  resistance,  R0  (b)  constant  s,  C  and  R0, 
and  various  values  of  the  voltage-dependent  component  of 
resistance,  r;  and  (c)  various  sweep  rates  for  selected  C,  /?0 
and  r  values. 

Two  cases  each  for  s  =  20  mV  s-1,  C  =  1  F  and  Rq  =  1, 
5,  10,  25  and  50Q  are  considered:  the  first  taking 
r  =  IQ  V^1  in  Eq.  (6)  and  the  second  r  =  5  Q  V^1.  Plots 
of  reduced  current  as  a  function  of  measured  potential  for 
r  =  1  Q  V-1,  together  with  that  for  r  =  ()QV  1 ,  for  various 
Rq  values  are  shown  in  Fig.  2a.  The  complementary  plot  of 
reduced  current  as  a  function  of  local  electrode  potential  is 
shown  in  Fig.  2b. 


Fig.  2.  Calculated  cyclic  voltammograms  expressed  as  reduced  response 
current  ( I/sC)  as  a  function  of  (a)  measured  potential  (Ve  +  I(Rq  +  rVe)) 
and  (b)  local  working  electrode  potential  (Ve)  calculated  using  Eqs.  (8)  and 
(9)  for  C  —  1  F,  V\  =  0  V,  V2  =  1  V,  j  =  20mVs-1.  Solid  curves 
r=lQV-1  and  symbols  r  =  0QV-1;  ('ft')  Indicates  direction  of 
increasing  Rq  (1,  5,  10  and  25  Q);  (— >-)  direction  of  voltammetric  sweep. 


Upon  examination  of  Fig.  2a,  it  is  clear  that,  for  a  small  r 
to  R0  ratio  (r/R0  <  1  V-1),  the  anodic  scan  remains  essen¬ 
tially  unchanged  from  that  for  the  r  OQV'  1  case, 
whereas  the  cathodic  reduced  current  response  is  initially 
more  sloping  (i.e.  resistive).  Toward  the  less  positive  poten¬ 
tial  limit  of  the  cathodic  scan,  the  reduced  current  response  is 
again  superimposable  upon  the  response  for  r  =  0  Q  V'  1 . 
Note  that  the  reduced  current  response  for  r  =  1  Q  V-1  does 
not  exhibit  the  same  symmetry  as  for  the  r  =  0  Q  V-1  case. 
From  Fig.  2b,  it  is  clear  that  the  r  component  further  limits 
the  upper  voltage  window  of  operation,  and  thus  reduces  the 
power  the  electrode  is  capable  of  delivering. 

If  r  is  comparable  to  or  larger  than  R0  (Fig.  3a  and  b),  then 
both  the  anodic  and  cathodic  reduced  current  responses  are 
modified.  In  this  case,  the  anodic  side  of  the  CV  has  the  same 
general  shape  as  for  the  case  when  r  =  0  Q  V-1,  but  the 
maximum  reduced  current,  I/sC~l,  achieved  is  <1.  Again, 
the  cathodic  response  at  more  positive  potentials,  is  more 
sloping  and  interestingly,  the  reduced  current  at  less  positive 
potentials  can  be  greater  than  I  (Fig.  3  for  r  =  5  Q  V  1  and 
Ro  =  1,  5,  10  Q).  For  these  cases,  the  plots  of  reduced 
current  as  a  function  of  local  or  measured  potentials  are 


(a)  V  m  /  V 


(b)  Ve  /  V 

Fig.  3.  Calculated  cyclic  voltammograms  expressed  as  reduced  response 
current  {II sC)  as  a  function  of  (a)  measured  potential  (Ve  +  I(Rq  +  rVe)) 
and  (b)  local  working  electrode  potential  ( Ve )  calculated  using  Eqs.  (8)  and 
(9)  for  C  —  1  F,  Vi=0V,  V2  =  1  V,  j  =  20mVs_1.  Solid  curves 
r  =  5QV-1  and  symbols  r  =  0QV-1;  (■ft)  indicates  direction  of 
increasing  R0  (1,  5,  10  and  25  Q);  (— >-)  direction  of  voltammetric  sweep. 
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Fig.  4.  (a)  Charge  and  (b)  energy  as  a  function  of  Rq,  calculated  from  the  CV’s  shown  in  Figs.  2  and  3:  (#)  r  —  1;  (■)  r  —  5;  ( — )  charge  cycle;  ( - ) 

discharge  cycle. 


highly  asymmetric,  indicating  that  the  local  electrode 
sweep  rate,  dVe/dt,  can,  under  certain  conditions,  be  larger 
than  s. 

In  Fig.  4a,  the  charge  capacity,  Q  in  (As)  calculated  from 
the  CV’s  shown  in  Figs.  2  and  3.  The  charge  efficiency 
determined  as  (Qoutput/Qmput)  x  100  was  calculated  to  be 
constant  and  approximately  100%  for  all  values  of  R0  and  r. 
The  energy  realisable  on  discharge  of  the  capacitor  elec¬ 
trode,  especially  for  cases  of  large  R0 ,  will,  however,  be 
much  less  than  that  involved  in  charging  the  same  electrode, 
thus  corresponding  to  serious  dissipative  losses.  This  is 
clearly  shown  in  Fig.  4b  in  which  the  energies  required 
for  charge  and  delivered  on  discharge  are  plotted  as  a 
function  of  R0.  The  energy  efficiency,  of  course,  decreases 
with  increasing  R0.  For  r/R(l  <  0.5  V  1 ,  the  absolute  max¬ 
imum  reduced  anodic  or  cathodic  current  densities  are  <  1 , 
and  the  CVs  become  more  symmetrical. 

In  Fig.  5a  and  b,  the  effect  of  increasing  the  ratio  of  r/R0 
on  the  reduced  current  as  a  function  of  reduced  and  local 
electrode  potentials,  respectively,  is  examined  for  a  given 
constant  sweep  rate  of  20  mV  s_1,  a  constant  capacitance  of 
1  F,  a  fixed  and  significant  value  of  Rq  =  10  Q  and  r/R0 
values  of  0, 0.1,  1,  2,  5  and  7.5  V-1.  In  Fig.  5c,  the  calculated 
charge  capacity  (i.e.  charge,  Q)  and  the  energy  required  for 
charge  and  that  delivered  on  discharge  for  the  CV’s  shown  in 
Fig.  5a  and  b. 


Again,  the  CV  response  is  not  symmetric  about  the 
II sC  =  0  line,  and  the  degree  of  asymmetry  increases  with 
increasing  r  value.  As  r  is  increased,  the  anodic  scan 
becomes  more  sloping  and  the  maximum  anodic  current 
achieved  decreases  as  does  the  total  charge.  The  cathodic 
scan  response  is  more  complicated.  Increasing  r,  initially 
causes  the  cathodic  response  to  be  more  sloping,  and  the 
maximum  reduced  current  to  be  >1,  further  increase  of  r 
results  in  an  even  more  sloping  CV  response,  and  reduced 
currents  <1.  A  plot  of  the  reduced  current  versus  local 
potential  (Fig.  5b)  at  the  working  electrode  illustrates  that 
increasing  r  reduces  the  potential  range  over  which  the 
working  electrode  actually  cycles. 

The  effect  of  reducing  the  potential  range  over  which  the 
working  electrode  cycles  is  shown  quantitatively  in  Fig.  5c 
in  which  electrode  charge  capacity  is  plotted  as  a  function  of 
r.  In  fact,  Q  is  a  linear  function  of  the  magnitude  of  the 
potential  range  over  which  the  working  electrode  actually 
cycles,  AVe,  as  shown  in  the  inset  of  Fig.  5c.  Although  the 
charge  capacity  decreases  with  r,  the  charge  efficiency  is 
independent  of  r  and  approximately  100%.  Both  the  energy 
required  for  charge  and  that  delivered  on  discharge  decrease 
with  increase  in  r  (Fig.  5d),  as  does  the  energy  efficiency. 

Finally,  the  effect  of  sweep  rate  on  CV  response  for  a 
system  defined  by  C  =  IF,  Rq  =  5  Q,  r  =  10  Q  V-1  and 
s  =  10,  20,  50,  100  and  200  mV  s-1  is  shown  in  Fig.  6.  At 
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(a)  Vm/V  (b)  VJW 


Fig.  5.  Calculated  cyclic  voltammograms  expressed  as  reduced  response  current  ( I/sC)  as  a  function  of  (a)  measured  potential  (Ve  +/(/?o  +  **Ve))  and  (b) 
local  working  electrode  potential  (Ve)  calculated  using  Eqs.  (8)  and  (9)  for  C  =  1  F,  V\  —  0  V,  V2  —  1  V,  s  =  20  mV  s-1,  Ro  —  10  Q,  symbols  r  =  0  Q  V-1. 
Solid  curves  r  —  1,  10,  20,  50  and  75  Q  V-1;  (It)  indicates  direction  of  increasing  r;  (— >-)  direction  of  voltammetric  sweep,  (c)  Charge  and  (d)  energy  as  a 
function  of  r  calculated  from  the  CV’s  shown  in  (a)  and  (b):  (♦)  charge  cycle;  (#)  discharge  cycle,  (c)  Charge  as  a  function  of  AVe  from  CV’s  shown  in  (a) 
and  (b). 


low  sweep  rate  (10,  20  and  50  mV  s-1)>  the  cathodic  current 
response  deviates  from  the  ideal  rectangular  response  much 
more  than  does  the  anodic  response.  Consider  first  the 
anodic  scan  at  low  sweep  rate;  at  low  potential,  the  local 
electrode  potential  is  approximately  the  same  as  the 
measured  potential,  and  the  electrolyte  resistance  is  simply 
R0.  At  higher  potential,  however,  the  rVc  component  of  the 
resistance  becomes  significant  resulting  in  a  limiting  value 
for  reduced  current  <1.  This  IrVe  component  then  results 
in  a  sloping  cathodic  current  response  which  causes  the 
cathodic  response  over  the  full  potential  range  to  be  far 
from  ideal. 

As  sweep  rate  increases,  IR  on  both  anodic  and  cathodic 
sweeps  increases  and  both  anodic  and  cathodic  responses 
deviate  from  the  rectangular  ideal  form,  (s  =  100  and 
200  mV  s-1)-  Further,  the  voltage  range  over  which  the  work¬ 
ing  electrode  operates  becomes  greatly  reduced  (Fig.  6b),  as 
does  the  charge  capacity  of  the  electrode  which  falls  from  1  C 
at  s  =  0  mV  s_1,  to  0.12  C  at  s  =  200  mV  s-1.  Again,  the 
charge  efficiency  independent  of  .v  and  100%. 

In  Fig.  6d,  the  plot  is  of  energy  to  charge  and  the  energy 
recovered  on  discharge  as  a  function  of  .v.  Both  quantities  are 
found  to  decrease  with  increase  in  s  (i.e.  decrease  with 
increase  in  charge/discharge  rate),  as  is  the  energy  efficiency 
because  increase  in  five  results  in  an  increase  in  both  the  IRq 
and  IrVe  losses. 


5.2.  Calculated  constant  current  charge/discharge 
behaviour 

Three  cases  are  considered  for  a  capacitative  electrode 
having  a  capacitance  of  1  F  operating  within  a  voltage 
window  of  1  V:  (a)  constant  I.  r  and  various  R0;  (b)  constant 
I,  R0  and  various  r;  and  (c)  constant  r,  R0  and  various 
constant  currents. 

5.2.1.  Case  (a) 

As  Rq  increases  the  “IR”  step  which  arises  at  the  begin¬ 
ning  of  charge  and  discharge  increases  as  seen  in  the  Vm 
versus  t  plots  of  Fig.  7a  for  second  and  subsequent  charge/ 
discharge  cycles  of  an  electrode  having  C  =  1 F, 
AV(measured)  =  1  V,  /  10  mA,  r  =  5  Q  V-1  and 

Rq  =  0,  1,  5,  10  and  25  Q.  As  a  result,  the  total  time  of 
charging  to  a  given  voltage  limit  (Vm  =  V2)  decreases  with 
increased  R(]  (therefore,  the  electrode  has  less  effective 
charge  capacity  within  a  particular  voltage  window).  Also, 
the  average  voltage  on  charge  is  greater  than  the  average  on 
discharge;  as  a  result,  power  consumption  during  charging 
will  be  greater  than  power  available  on  discharge  (i.e.  watt- 
hour  efficiency  will  be  low,  and  decrease  with  increasing 
Rq).  a  significant  practical  conclusion.  In  Fig.  7b,  the  local 
electrode  potential,  Ve,  as  a  function  of  time  during  charge/ 
discharge  between  measured  potential  limits  of  0  and  1  V. 
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Fig.  6.  Calculated  cyclic  voltammograms  expressed  as  reduced  response 
current  ( I/sC)  as  a  function  of  (a  measured  potential  ( Ve  +  I(R0  +  rVe) ) 
and  (b)  local  working  electrode  potential  (Ve)  calculated  using  Eqs.  (8)  and 
(9)  for  C=1F,  V,=0V,  V2  =  1  V,  R0  =  10 r=10nv_1.  (f) 
Indicates  direction  of  increasing  s;  (— >-)  direction  of  voltammetric 
sweep,  (c)  Charge  and  (d)  energy  calculated  from  the  CV’s  shown  in  (a) 
and  (b):  (♦)  charge  cycle;  (#)  discharge  cycle. 


Clearly,  the  capacitative  electrode  is  cycling  within  a 
restricted  potential  range,  which  decreases  with  increasing 
R0.  Note  also  that  increasing  If  ,  both  decreases  the  anodic 
potential  limit  and  increases  the  cathodic  limit,  thus  redu¬ 
cing  the  cycling  window  at  both  extremes  of  potential. 

5.2.2.  Case  (b) 

Fig.  8a  and  b  plotted  the  Vm  and  Vc  versus  1,  respectively, 
for  an  electrode  defined  by  C  =  1  F,  Admeasured)  =  1  V, 
/=  10  mA,  Ro  =  10  and  rofO,  10, 20,  35, 50  and  75  QV1. 
The  IR  drop  on  switching  from  discharge  to  recharge  is 
determined  primarily  by  Ro  (except  in  the  case  of  very  large 
r/R0).  while  that  for  switching  from  charge  to  discharge  is 
due  to  the  combination  of  r  and  R0  terms.  Interestingly, 
however,  the  working  electrode  reversal  potential  on  the 
second  and  subsequent  charge  is  determined  by  the  magni¬ 
tude  of  r.  Again,  increasing  r  will  lower  both  watt-hour 
efficiency  (due  to  higher  charging  voltage  as  compared  with 
discharge  voltage)  and  charge  storage  capacity  (as  r 
increases  the  local  working  electrode  cycling  range 
decreases). 


Fig.  7.  Calculated  voltage  profile  for  second  constant  current  charge/ 
discharge  cycle  as  a  function  of  time  using  Eq.  (11)  for  C  =  1  F,  Vi  =  0  V, 
V2  =  1  V,  /  =  1 0  mA,  r  =  5  Q  V  1 :  (a)  measured  potential 

(Ve  +  /(Ro  +  rVe))  and  (b)  local  working  electrode  potential  (V'e).  (f[) 
Indicates  direction  of  decreasing  R0  (25.  10.  5,  1  and  0  Q). 


5.2.3.  Case  (c) 

The  rate  effects  for  a  system  having  C  =  1  F,  AV 
(measured)  =  1  V,  Rq  =  5  Q  and  r  =  10  V-1  and 

I  =  10,  20,  30  and  40  mA  are  shown  in  Fig.  9a  and  b  for 
measured  and  local  electrode  potentials,  respectively,  as  a 
function  of  charge  (mAh)  for  second  and  subsequent  charge/ 
discharges.  The  effective  electrode  charge  capacity  (Ah)  and 
effective  capacitance  (F)  are  both  rate  dependent.  Capacity 
(Ah)  decreases  with  increasing  I  (and  consequently  increas¬ 
ing  IR).  As  the  current  increases  the  effective  capacitance 
increases,  but  over  a  progressively  decreasing  voltage  win¬ 
dow.  The  charging  efficiency  (Ah)  remains  relatively  con¬ 
stant  with  increasing  I  (>90%).  Similar  to  the  results  shown 
in  Figs.  7  and  8,  increasing  rate  results  in  a  restriction  of  the 
local  potential  cycling  range.  Similar  trends  are  shown  in  the 
three  Figs.  7-9  where  it  is  seen  that  each  of  increasing  R0,  r 
and  I  results  in  an  increase  in  the  “IR”  loss  component. 

5.3.  Comparison  with  experimental  results  for  specially 
prepared  non-aqueous  carbon  electrode  supercapacitors 

5.3.1.  Cyclic  voltammetry  (CV) 

The  experimentally  derived  capacitance  response  voltam¬ 
mograms  for  specially  prepared  [18]  capacitor  devices 
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(a)  time  /  s 


Fig.  8.  Calculated  voltage  profile  for  second  constant  current  charge/ 
discharge  cycle  as  a  function  of  time  using  Eq.  (11)  for  C  =  1  F,  Vi  =  0  V, 
Wn  =  1  V.  /  =  10  mA,  Rt]  =  10  Q:  (a)  measured  potential 

(Ve  +  I(Ro  +  i'Ve))  and  (b)  local  working  electrode  potential  (14).  (if) 
Indicates  direction  of  decreasing  r  =  0.  10,  20,  35,  50  and  75  Q  V  . 

having  1 .0  and  0.4  M  TEATFB  electrolyte  are  shown  in 
Fig.  10a  and  b,  respectively,  as  a  function  of  sweep-rate. 

In  Fig.  10a,  for  the  1.0  M  electrolyte,  the  capacitance 
behaviour  approximates  to  ideal  behaviour  at  the  lowest 
sweep-rate  of  1  mV  s  1 ,  i.e.  a  voltammogram  having  an 
almost  rectangular  form  and  mirror-image  symmetry  of  the 
current  responses  about  the  zero-current  line  is  observed. 
With  increasing  s  the  voltammograms  become  distorted 
from  the  rectangular  form  expected  for  an  ideal  capacitor 
having  constant  capacitance,  similar  to  the  calculated  effects 
shown  in  Fig.  1. 

The  CV  response  of  the  module  made  with  0.4  M  elec¬ 
trolyte,  shown  in  Fig.  10b,  for  s  =  5,  10,  20,  30  and 
40  mV  s“  ,  exhibits  behaviour  similar  to  that  calculated 
for  the  significant  resistance  models  of  this  study  (cf. 
Fig.  2).  Here  the  capacitance  response  as  a  function  of 
voltage  is  nearly  ideal  at  low  sweep  rate  (5  mV  s-1),  while 
at  higher  sweep  rates,  becomes  distorted.  At  high  s,  unlike 
the  ideal  behaviour,  the  realised  capacitance  is  both  varying 
with  module  voltage  and  increasing  with  decreasing  s. 

5.3.2.  Constant-current  charge/discharge  behaviour 

The  observed  constant-current  charge/discharge  beha¬ 
viour  is  plotted  in  Fig.  11a  and  b  for  the  respective  capacitor 


Fig.  9.  Calculated  voltage  profile  for  second  constant  current  charge/ 
discharge  cycle  as  a  function  of  time  using  Eq.  (11)  for  C  =  1  F,  Vi  =  0  V, 
V2  =  1  V,  So  =  5  fi,  r  =  10  fi  V  1 :  (a)  measured  potential 
(Ve  +  1(A)  +  rVc))  and  (b)  local  working  electrode  potential  (14).  (If) 
Indicates  direction  of  decreasing  I  =  10,  20,  30  and  40  mA. 

modules  discussed  in  Section  5.3.1,  having  1.0  and  0.4  M 
electrolyte,  respectively,  as  a  function  of  calculated  cumu¬ 
lative  charge,  Q ,  in  As  at  constant  current  charge/discharge 
rates  of  1.33,  1.0,  0.67,  0.33  Ag_1  of  carbon. 

Note,  even  for  the  1.0  M  module  shown  in  Fig.  11a,  that 
there  is  a  small,  but  significant  sharp  drop  in  voltage  on 
changing  from  charging  to  discharging  current  and  a  sharp 
increase  at  the  start  of  the  charging  process.  This  is  not  due 
entirely  to  a  direct  ESR  effect  but  is  also  a  characteristic 
of  redistribution  of  charge  within  the  pores.  The  same 
characteristic  curve  was  demonstrated  in  [22]  found  for  a 
hardware,  four-element  model  RC  network,  made  with 
relatively  small  resistive  elements.  Probe  measurements 
enabled  redistribution  of  charge  at  each  of  the  capacitor 
elements  in  this  hardware  RC  ladder  network  to  be  directly 
observed  experimentally  [22]. 

Fig.  lib  shows  similar  plots  for  the  capacitors  employing 
the  smaller  concentration  of  electrolyte,  0.4  M.  Here  the  IR 
drops  are  significantly  larger  than  those  in  Fig.  11a,  and  even 
at  the  lowest  rate  of  charge/discharge  of  0.25  A,  the  charge 
capacity  of  the  0.4  M  electrolyte  module  is  significantly 
reduced  from  that  of  the  1.0  M  module.  This  further  illus¬ 
trates  that  in  the  low  concentration  electrolyte  modules  the 
local  electrode  potential  range  is  much  smaller  than  the 
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Fig.  10.  (a)  Cyclic  voltammograms,  showing  the  current  responses 
expressed  as  capacitance,  for  the  1 .0  M  TEATFB  electrolyte  module  as  a 
function  of  voltage  over  the  voltage  range  0-3  Vat  1,  5.  10  and  20  mV  s  1 ; 
arrow  indicates  direction  of  increasing  sweep  rate;  (b)  cyclic  voltammo¬ 
grams,  showing  the  current  responses  expressed  as  capacitance,  for  the 
0.4  M  TEATFB  electrolyte  module  as  a  function  of  voltage  over  the 
voltage  range  0-3  V  at  5,  10,  20,  30  and  40  mV  s  1 .  ("ft")  Indicates  direction 
of  increasing  sweep  rate. 

measured  range.  Thus,  IR  effects  can  result  in  reduced 
energy-  and  power-density,  and  therefore,  performance 
characteristics  of  practical  devices.  A  rate  increase  of 
0.33-1.33  A  g-*1  of  carbon  results  in  a  reduction  of  the 
charge  capacity  from  0.99  to  0.72  As  (30%  reduction)  for 
the  1.0  M  electrolyte  device,  and  from  0.61  to  0.21  As  (66% 
reduction)  for  the  0.4  M  electrolyte  device.  This  results  from 
a  restriction  in  the  local  electrode  potential  range  as  pre¬ 
viously  discussed,  and  leads  to  the  practically  significant 
reduction  in  energy  storage  capacity  and  energy-density. 

These  features  remain  largely  unrevealed  in  constant- 
current  charge/discharge  curves  for  well  engineered,  prac¬ 
tical  capacitors  having  an  optimised  low  resistance  of  elec¬ 
trolyte,  especially  in  the  cases  of  aqueous  electrolyte 
embodiments. 


6.  Conclusions 

Increasing  s ,  r,  or  R0  results  in  a  decrease  in  the  overall 
charge  storage  capacity  of  the  electrode.  The  ampere- 


Fig.  1 1 .  The  voltage  response  to  constant  current  charging/discharging 
currents  of  1.33,  1.0,  0.67  and  0.33  A  g_1  of  carbon,  as  a  function  of 
cumulative  calculated  charge  (As)  for  a  redox  engineering  3  V  AA 
electrochemical  capacitor:  (a)  electrolyte  concentration  1.0  M,  and  (b) 
electrolyte  concentration  0.4  M.  (1))  Indicates  the  direction  of  decreasing 
current  density. 


hour  efficiency  remains,  however,  constant,  at  100%  for 
all  s ,  r  and  R0  considered  (i.e.  ampere-hour  input  =  ampere- 
hour  output).  The  energy  efficiency  (watt-hour  output/ 
watt-hour  input)  was  found  to  decrease  with  increasing  s, 
r  or  Rq.  Any  variation  in  s,  r  or  R0  which  increases  the 
asymmetry  of  the  CV  response,  results  in  a  decrease  in  the 
available  energy,  and  thus  a  decrease  in  energy-density 
an  operational  energy  efficiency.  The  modelling  treatment 
provides  a  good  computational  basis  for  interpreting 
the  experimental  behaviour  of  selected  capacitor  test 
modules. 
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